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Summary. — We propose a solution of the cusp problem in framework of the 
standard ACDM cosmology. To do this we describe the linear and nonlinear periods 
of halo formation by the entropy function of dark matter particles. This approach 
allows us to take into account together the impact of both the processes of nonlinear 
relaxation of compressed matter and the small scale initial velocity perturbations in 
collapsed halos. We show that such random velocities lead to the random variations 
of the density profile of relaxed halos. As a rule, they suppress the formation of 
cusp-like halos and favor the creation of core-like ones. This approach allows us to 
reproduce observed rotation curves, to explain their random scatter and deviations 
from simulated ones. 
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1. Introduction 

The cusp problem comes from numerical simulations predicted the formation of di- 
vergent density profiles inside virialized DM halos. However such a behavior of density 
profiles is not seen in observations. This contradiction is known as "the cusp problem". 

Mathematically this means that in the central part of gravitationally-bounded halos 
we get the powerdaw density profile 



(1) 



p{r) oc r 



with a > 1 in numerical simulations (e.g. [Tj) and a < 1 in observed galaxies (e.g. [21 [3]). 
In the first case we deal with cusp, in the second one - with core. In addition, the available 
observational data show a broad scatter of a (see [H [3]) and sometimes a > 1 are 
observed (for example, in rich clusters of galaxies [4]). Majority of the observed rotation 
curves of galaxies [5] are well fitted by the Burkert's function with a — (see [5]) which 
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differs significantly from the simulated curves fitted by the Navarro-Frenk- White (NFW) 
function with a = 1, ([3 [3]). 

This situation has created the problem discussed during many years as an important 
problem of the standard ACDM cosmology. More of that, there is a common myth that 
the relaxation of collisionless DM particles always produces density cusps. 

2. Physical nature of difference between cusps and cores 

To clarify the difference between core and cusp we calculate pressure profile for DM 
halo. For power-law density profile (JXJ) we get for the mass of halo 

/•OO 

(2) M = Mir) = \ p(r)r 2 dr oc r (3 " a) , 

Jo 

and for the pressure, p(r), within an equilibrium halo 

ldp_ GM{r) (i-g) 
p dr r z 



(3) p{r)=C 1 +C 2 r 2( - 1 - a \ 

where C\ and C 2 are some constants. So, if a < 1 then in the central part of halo the 
pressure is finite and we have the core case. If a > 1 then the pressure is infinite what 
corresponds to the case of cusp. Evidently, a — 1 is a critical value of a corresponding 
to the logarithmic cusp formation. 

The halo formation includes both processes of the reversible matter compression and 
the irreversible relaxation. To discriminate between them we will consider the entropy 
function, F(r). For the nonrelativistic DM particles with the isotropic pressure it can be 
introduced as follows: 

(4) p(r) = p(v 2 ) =nT = F(r)n 5/3 (r) , 

where v is one-dimensional random velocity, p, n, & T are the density, concentration and 
temperature of DM particles. F can be presented as a function of current radius of halo, 
r, but it is more convenient to consider F as a function of current mass, i.e. mass inside 
of r. We call F(M) by entropy mass function: 

(5) F(M) oc CiM 131 + C 2 M f)2 oc M 13 , 



(6) h = vo a v ^ 2 = ^rA' 

6(6 — a) 3(3 — a) 

For the critical value of a = 1 we get j3 cr = (3± = fa = 5/6. 

The entropy mass function integrates the impact of irreversible processes during all 
DM evolution and determines the profile of relaxed halo. As is seen from © for a < 3 
we have j3\ > 0, fa > and, so, F[M) — > for M — ► 0. However, the fraction of low 
entropy DM particles in the central regions of halo in the cusp case is larger than that 
in the case of core. 
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3. Idea and method 

The discussions of the halo formation are usually concentrated around the processes 
of DM relaxation which are main sources of entropy. At the same time it is usually 
accepted that before matter relaxation the DM entropy is negligible. However, in the 
CDM model there are random correlated velocities down to very small scales determined 
by the mass of DM particles. In a course of the matter compression they are partly 
transformed to the temperature of particles and, so, increase the entropy mass function 
at small M. This effect can be estimated analytically. Further on we will calculate the 
joint entropy in relaxed DM halos summing the initial entropy given by small scale initial 
perturbations and an entropy generated during non-linear relaxation of collapsed matter. 

Our further analysis includes the following steps: 

1. Instead of density-radius relations we use the entropy-mass profiles. 

2. The mean initial entropy profile, Fb(M), is estimated as the entropy mass function of 
linear velocity perturbations. Such identification assumes the almost adiabatic matter 
compression and, so, partly underestimates the entropy mass function i*&(M) at small 
M. 

3. The entropy generated by the process of relaxation of compressed matter, F r (M), is 
taken from results of N-body simulations (e.g. NFW profiles). 

4. We combine the initial and generated entropies by a simple way to set of joint entropy 
profiles and reconstruct the circular velocities and resulting density profiles. 

4. Initial entropy mass function 

At the beginning we need to consider initial field of density perturbations and its 
statistical properties. 

Let us characterize the initial conditions in dark matter by three gauge-invariant 
variables - displacement of a matter point from unperturbed position (the deviation 
from Hubble flows), S(x), the full velocity of matter, V(z, x), and comoving density 
perturbation, 6(z,x) = p(z,x)/p(z) — 1: 

x-g(z)S(x) , 
(76) V{z,x)='r = H(r + g'S), 

8{z 1 x) = g(z) ■ S(x) , 5(x) = div S = dSi/dxi , 

where the dot and prime mean derivative with respect to the physical time and the 
function argument, f(z, x) and x are the Eulerian (shear- free) and Lagrangian (comoving) 
coordinates of a matter point. 

In these terms perturbations related with the collapsed halo as whole and conditional 
perturbations within the halo can be separated as follows: 

(8a) S(x) = S R (x) + S*(x) 

(86) 5(x) = S R {x) + S*(x) 

where quantities with index R are related to a protohalo with a linear size R, and 
quantities with index * are related to conditional perturbations which determine the 



(7a) 



f(z, x) 



1 4- z 
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Table I. - Power index of the function F oc M 13 for halos with M < 1O"M . 



n = log(M/M Q ) 


10 


7 


5 


< 5 


p = d\nF{M)/d\nM 


0.3 


0.5 


0.57 


0.67 



corresponding temperature and entropy of compressed matter, 

(9) (&) = (ft,) = (S 2 ) = al(R) = J P(k) [1 - W{kR)fdk. 

where W{x) is a usual spherical window function. 

After some analytical and numerical calculations we obtain for the power index /? of 
initial entropy mass function F(M) for M < 10™M Q (see Table H}: 

It is easy to see that in all cases is less than its critical value j3 cr = 5/6. This means 
that initial entropy dominate in central regions of halo what stimulates the formation of 
cores and suppresses the cusp formation. 

5. Entropy generation during the matter relaxation 

The spherical collapse with the violent relaxation of compressed matter have been 
investigated in many papers. Thus, authors of [9] started from initial conditions with no 
initial velocities and entropy 

(10) SM(r) = 1 - M(r)/{M(r)) cx {M(r))- £ , 

(11) v t = 0, Sp = p m (r) - (p in ) ocr 3e (M(r))cxr 3 , 

and found that the density profile of relaxed matter, p(r), is approximated at small r by 
power law: 

(12) P{r)<xr-", a=| 9£/{l \ 3eh J | 2/3 ' 

The approach was extended in [101 IH] where nonradial trajectories of DM particles 
were also considered. The density profiles for relaxed halos was found to be similar to 
previous one. 

This problem have also considered in [T2J [T3j for the initial conditions 

(13) «(r) = 0, P {x) = p (l - r 2 /r 2 ) . 

Near the center of cloud the power law density profile was found with 

(14) p(r)(xr- a , a ~ 1.7 -1.9. 

Anywhere initial velocities and entropy was accepted to be equal to zero. 



A SOLUTION OF THE CUSP PROBLEM IN RELAXED HALOS OF DARK MATTER 



5 



As it is seen from these relations, the model of violent relaxation leads to the entropy 
distribution 

(15) M(r)ocr, F g (r) oc r 4/3 oc M 4/3 , fi = 4/3 > fi cr , 

that implies the formation of cusp. However, in the central regions of halos the entropy 
generated by this process becomes negligible as compared with the initial entropy. Thus, 
the influence of former one restricts the central density of halos and prevents formation 
of cusp-like density profile. 

The same conclusion is valid also for the simulated NFS - density profile with a — 
a cr = 1, fi = (3 cr = 5/6 which represents the more general model of hierarchical halo 
formation. The drop of these indices as compared with previous ones indicates that in 
the case the initial entropy is partly allowed for but owing to the small scale cutoff of 
the simulated power spectra its impact is underestimated as compared with estimates in 
Table 1. 

6. Rotation curves caused by the joint entropy 

The shape of rotation curves is determined by sum of initial and generated entropy. 
To clarify influence both of them we model the joint entropy mass function as follows: 

(16) F{M) = Jf[+f], F b = KF. 

where the random factor k, < n < 1, measures the relative contribution of the initial 
entropy for different halos. 

Now we can calculate rotation curves and compare them with observed and simu- 
lated ones. The results are presented in Fig. 1 where the circular velocities, v c , are 
plotted versus radius for six models of hierarchal (fi g = 5/6, left plot) and six mod- 
els of violent (fi g = 4/3, right plot) relaxation processes for initial conditions with 
fib — 0.333,0.567,0.667 and for k < 1 and k ~ 1. Both the velocities and radius 
are normalized for their values in point of velocity maximum, v max — v c {r max ) > v c (r). 
The curves are compared with the NFW and Burkert ones. 

As is seen from this Figure, the rotation curves for our models with fib — 1/3 and 
< n s < 1 close the gap between the Burkert and NFW rotation curves. So, we can 
expect that properties of observed curves can be successfully reproduced by our models 
with suitable parameters k and fib- The scatter of both initial and generated entropies 
can provides required variations of these parameters and the shape of rotation curves 
what, in turn, explains variations of observed rotation curves. 

Let us note also that for low mass halos with [3b — 0.567 and especially fib — 2/3 
all rotation curves with < n < 1 are concentrated nearby the NFW profile. This fact 
indicates that probably for dwarf galaxies with M < 10 8 M Q the observed rotation curves 
can be similar to expected ones for the model of hierarchical clustering with the NFW 
profile. 

7. Conclusions 



We conclude, that 
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Fig. 1. - Normalized rotation curves are plotted for the models of hierarchical (left side of the 
Figure) and violent (right side) relaxation for k 3 <C 1 and k s w 1 (dashed and solid lines) . NFW 
and Burkert fits are plotted by stars and dots, respectively. 



1. The initial entropy can prevent the cusp formation for halos with DM masses 
10 8 — 10 12 M Q . For smaller and larger galaxies and for clusters of galaxies the impact of 
the initial entropy is attenuated. 

2. The impact of the initial entropy allows to reproduce the observed rotation curves 
and many helps to solve the so called "cusp problem" . 

* * * 
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